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Abstract —The phasor measurement unit (PMU) placement 
problem is revisited by taking into account a stronger charac¬ 
terization of the electrical connectedness between various buses 
in the grid. To facilitate this study, the placement problem 
is approached from the perspective of the electrical structure 
which, unlike previous work on PMU placement, accounts for 
the sensitivity between power injections and nodal phase angle 
differences between various buses in the power network. The 
problem is formulated as a binary integer program with the 
objective to minimize the number of PMUs for complete network 
observability in the absence of zero injection measurements. The 
implication of the proposed approach on static state estimation 
and fault detection algorithms incorporating PMU measurements 
is analyzed. Results show a significant Improvement in the 
performance of estimation and detection schemes by employing 
the electrical structure-based PMU placement compared to its 
topological counterpart. In light of recent advances in the 
electrical structure of the grid, our study provides a more realistic 
perspective of PMU placement in the electric power grid. 

Index Terms —PMU placement, topology, electrical structure, 
state estimation, fault detection. 

I. Introduction 

In the classic setting, the problem of placing PMUs within 
an electric power system is divided into two parts: 1) obtaining 
the optimal or minimum number of PMUs; and 2) hnding 
the optimal locations to install these PMUs on the grid. In 
both stages, the optimizations are typically done with the goal 
to have either complete or incomplete network observability 
either in the presence/absence of zero injection measurements 
m . Zero injection measurements are present when the power 
system has nodes without generation or load. The placement 
problem is formulated as an integer linear program, comprising 
the adjacency (or connectivity) matrix of the grid, to obtain 
the minimal set of PMUs and their optimal locations. 

Various strategies have been investigated to address the 
PMU placement problem. For instance, in 12], the placement 
was provided by a spanning measurement subgraph, and the 
minimal PMU set was obtained using a dual search algorithm 
incorporating modified bisecting simulated-annealing-based 
searches for complete network observability. In fS), a Pareto- 
optimal solution was obtained using non-dominated sorting 
genetic algorithm, which turned out to be nonlinear in the 
presence of power injection measurements. In Q, the problem 
was formulated as an integer program to include conventional 
power flow and injection measurements in addition to PMU 
measurements for maximum network observability. In Q, 
the objective was to achieve bad data detection during state 
estimation using optimal PMU placement. In l6|, I?], optimal 
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placements were obtained by formulating the problem as a 
simple integer linear program to achieve standard objectives. 

In ISl, the optimal placement was linked with power system 
dynamic state estimation. A method for evaluating a specific 
placement configuration, when there were multiple solutions, 
was also proposed. In a, the problem was solved using 
an exhaustive search-based method with complete network 
observability, and the state estimation implemented on such a 
placement was shown to be linear. A convex relaxation method 
was employed in ca to optimize the placement based on 
estimation-theoretic criteria under the Bayesian framework. In 
na, an information-theoretic measure, namely, mutual infor¬ 
mation (MI) was employed to address the PMU placement 
problem, where the objective was to maximize the MI between 
PMU measurements and power system states to obtain highly 
“informative” PMU configurations. 

Studies in the aforementioned work were based on the 
topology or node connectivity (degree distribution) of the 
grid. However, advances in complex networks have questioned 
the veracity of characterization of the power grid using its 
topological structure ifT^ . Subsequently, several demerits of 
the topology-based analysis were uncovered. For instance: 

(1) Studies showed that, for many real-world power grids, 
characterizing the network structure using topology alone 
was suboptimal, and had implications on node synchro¬ 
nization and performance of the network ifTSll - ifTSll . 

(2) It was reported that grids in different geographical re¬ 
gions had different degree distributions, leading to varied 
topological structures. Furthermore, it was also shown 
that different model-based analyses of the same grid had 
resulted in different topological structures flhll . ifTTll . 

(3) It was shown that the topological structure did not com¬ 
prehensively capture the underlying physical laws (Ohm’s 
and Kirchoff’s) that govern the electrical connections or 
flows between network components lIT^ . 

(4) In the context of PMU placement, the topology-based 
approach provides a weak characterization of the electrical 
influence between various components in the grid, since 
it does not take into account the effect of bus voltage 
magnitudes and phasor angles. Also, owing to the the 
time-varying (albeit slowly) nature of the networks’ con¬ 
nectivity, an optimal PMU placement (minimum number 
and locations) apposite to a particular instant in time might 
not hold for a different time instant. 

To overcome the aforementioned drawbacks of the topological 
structure of the grid, authors of ifTSll studied the electrical 
structure, which accounted for the sensitivity between power 
injections and nodal phase angle differences between buses in 
the grid, to provide a more comprehensive characterization of 
the connectedness between various network components. 


Motivated by these observations, the PMU placement prob¬ 
lem is revisited in this paper from the perspective of the electri¬ 
cal structure of the grid. To the best of the our knowledge, this 
is the first instance of the placement problem being addressed 
from the perspective of the electrical structure, which provides 
a much stronger characterization of the electrical influence 
between various buses in the gird. As the results demonstrate, 
a larger number of PMUs are required when compared with 
the number required when using the topological structure to 
meet the objective of complete network observability in the 
absence of zero injection measurements. However, in light 
of the aforementioned arguments which favor the electrical 
structure of the grid, our results provide a more realistic 
perspective of PMU placement. The studies in this paper are 
topical especially with the proposition by the electric utilities 
to install large numbers of PMUs on the grid. 

The implication of the proposed approach on the perfor¬ 
mance of (i) static state estimation and (ii) fault detection 
algorithm to detect changes in the susceptance parameters of 
the grid is analyzed. The detection problem is formulated using 
a linear errors-in-variables (EIV) model, and a generalized 
likelihood ratio test (GLRT) based on the total maximum 
likelihood (TML) methodology is presented. Results show a 
significant improvement in the performance of both estimation 
and detection schemes by employing the electrical structure- 
based PMU placement compared to its topological counterpart. 

In Section |II] a review of the electrical structure of the grid 
is provided to build the necessary framework to be used in 
the rest of the paper. In Section |III1 the placement problem is 
formulated and the results are tabulated for several IEEE bus 
systems. State estimation and fault detection frameworks with 
PMU measurements are presented in Section |IV] Section m 
comprises experimental results and related discussion. We 
conclude the paper in Section |VT] 

II. Electrical Structure of the Power Network 

In this section, the electrical structure of the grid from the 
perspective of complex networks is reviewed. The concept of 
resistance distance is introduced, and the procedure to derive 
the binary connectivity matrix of the power grid is presented 
(see ifTSl Section III] for further details). 

The sensitivity between power injections and nodal phase 
angles differences can be utilized to characterize the electrical 
influence between network components. The electrical struc¬ 
ture of the power network can then be understood by measur¬ 
ing the amount of electrical influence that one component has 
on another in the network. The measurement of this electrical 
influence necessitates a metric system. Mathematically, this 
can be accomplished by first deriving the sensitivity matrix, 
which can be obtained by well-established methods. The 
complement of the sensitivity matrix is called the distance 
matrix, whose entries quantify the electrical influence that 
each component has on the other - zero value indicates that 
two components are perfectly connected, while a large number 
indicates that the corresponding components have negligible 
electrical influence on each other. This electrical distance was 
proved to be a formal distance metric (see iTTSl Appendix]), 
and was employed to address various problems in the grid. 


Another method to measure the electrical influence between 
network components is to derive the resistance distance on, 
which is the effective resistance between points in a network 
of resistors. Consider a network with N nodes, described by 
the conductance matrix G. Let Vj and gij denote the voltage 
magnitude at node j and the conductance between nodes i 
and j, respectively. The current injection at node i is given 
by Ii = ^ ^ Laplacian matrix to the 

network, provided there are no connections to the ground, i.e., 
if G has rank — 1. The singularity of G can be overcome 
by letting a reference node r have Vr = 0. The conductance 
matrix associated with the remaining — 1 nodes is full- 
rank, and thus 14 = G^^Ik,k ^ r, where is the sub¬ 
matrix of G associated with the non-reference nodes, and 14 
and Jfc denote the vectors of voltage magnitudes and current 
injections, respectively, at the non-reference nodes k^r. 

Let the diagonal elements of G^ be denoted V/c, 
indicating the change in voltage due to current injection at 
node k which is grounded at node r. The voltage difference 
between a pair of nodes i ^ j ^ r, is e{i,j) = 

+ gj^ - - gj^, indicating the change in voltage due to 

injection of 1 Ampere of current at node i which is withdrawn 
at node j. e{i,j) is called the resistance distance between 
nodes i and j, and describes the sensitivity between current 
injections and voltage differences. Letting T = diag(Gj4^) 
and 1 denote the vector of all ones, we have 


Ekk 


(1) 

Ej-k 

= r^, 

(2) 

^kr 

= r. 

(3) 


([Til - (O holds for Vfc ^ r. The resistance distance matrix E, 
thus defined, possesses the properties of a metric space El. 

To derive the sensitivities between power injections and 
phase angles, we start with the upper triangular part of the 
Jacobian matrix obtained from the power flow analysis, for 
the distance matrix to be real-valued; 


AP = 


'dP' 


A6»-f 


■ dP ■ 

W\. 


A|U|. 


(4) 


The matrix [^] will be used to form the distance matrix 
by assuming the voltages at the nodes to be held constant, 
i.e., A|U| = 0. It was observed that [^] possesses most 
of the properties of a Laplacian matrix (see El). By letting 
G — [^], the resulting distance matrix E measures the 
incremental change in phase angle difference between two 
nodes i and j, {9i — 9j), given an incremental average 
power transaction between those nodes, assuming the voltage 
magnitudes are held constant. It was proved in El Appendix] 
that E, thus defined, satisfies the properties of a distance 
matrix, as long as all series branch reactances are nonnegative. 

Eor a power grid with N buses, the distance matrix E 
translates into an undirected graph with N{N — 1) weighted 
branches. In order to describe a grid with an undirected 
network without weights, one has to retain the N buses, but 
replace the K branches with K smallest entries in the upper 
or lower triangular part of E. This results in a graph of size 
{N,M} with edges representing electrical connectivity rather 








than direct physical connections. The connectivity matrix C 
of this graph is obtained by setting a threshold, r, adjusted to 
produce exactly K branches in the network; 

^ fcy = 1, Ve(i,j) < r, 

■ jcij = 0, Ve(i, j) > T 

To obtain the threshold t in (|5]l, we first consider the upper 
triangular part of the matrix E, and sort the elements in 
descending (or ascending) order. We then pick a number of 
elements equal to the number of branches in the given power 
network. For instance, in the IEEE-14 bus network having 
20 branches, we pick the top 20 sorted elements from the 
upper triangular part of the matrix E. In this paper, the binary 
connectivity matrix C is derived for several standard IEEE 
test bus systems for use in the PMU placement problem. 

III. Problem formulation and main results 

The PMU placement problem is formulated and addressed 
from two different perspectives: one based on the topological 
structure of the grid and the other based on its electrical 
structure. The problem is setup as an integer linear program 
to achieve complete network observability in the absence of 
zero injection measurements. The extension to the case of 
incomplete network observability (see HI) is straightforward 
and, therefore, is not considered in this paper. 

Consider a power network with N buses and K branches. 
Let i and j to denote the bus indices, i,j G {1,. .., TV}, and 
let C denote the N x N binary connectivity matrix. Let d be 
an iV X 1 binary decision variable vector defined by 

1, if a PMU is installed at bus i, 

(6) 

0, otherwise, 

where i = 1,... , N, and b a unit vector of dimensions x 1. 
The PMU placement problem for complete network observ¬ 
ability and without conventional measurements is formulated 
as follows ©, Q: 

N 

min di 

i=l 

such that Cd ^ b. (7) 

Eor PMU placement based on the topological structure of 
the power network, we consider the existing approach of 
deriving the binary connectivity matrix C directly from the bus 
admittance matrix. The entries of the bus admittance matrix 
are transformed into binary form, and used in the problem 
setup ([7|i. The entries of C are given by 

{ dj = 1, if i = j, 

Cij = 1, if i and j are connected, (8) 

Cij = 0, if i and j are not connected. 

The entries dj characterize the electrical influence between 
the buses i and j based on the degree distribution of the grid. 
Eor the electrical structure-based approach, the connectivity 
matrix C is given by 0. 

The bus and branch data required to derive the bus ad¬ 
mittance and power-flow Jacobian matrices were obtained 


using MATPOWER Il20ll . Eor the topology-based approach, the 
connectivity matrix is computed using 0, which is employed 
in 0 to obtain the placement soution. Similarly, for the 
electrical structure-based approach, the connectivity matrix 
is obtained using © as described in Section El Eor both 
approaches, the criterion for the optimization problem was to 
achieve complete network observability, without conventional 
measurements. The binary integer programming tool of Matlab 
was used to solve the problem defined by ©. 


IEEE bus system 

Topological structure 

Electrical structure 

9 

3 

4 

14 

4 

7 

30 

10 

20 

39 

13 

24 

57 

17 

39 

118 

32 

98 

162 

43 

131 


TABLE I; Minimum number of PMUs based on topological 
and electrical structures for IEEE test bus systems 

The main result of the paper is tabulated in Table 11 It 
can be seen that for each test bus system, the minimum 
number of PMUs obtained using the electrical structure is 
higher than that obtained using the topological structure, to 
meet the objective of complete network observability in the 
absence of zero injection measurements. Though this result 
is unfavorable from an economic standpoint, the arguments 
presented in Section H clearly support the employment of 
more PMUs for optimal system operation. Eurthermore, the 
optimal locations of the PMUs are also different for these 
two approaches. Eor instance, for the 14 bus system, the 
topology-based approach suggests to place the 4 PMUs on 
buses numbered 2, 6, 7 and 9, which is in agreement with I©, 
Q. Whereas, the electrical structure-based placement scheme 
suggests to place the 7 PMUs on buses numbered 1, 3, 8, 11, 
12, 13 and 14. A similar behavior was noticed for the other bus 
networks, however, in the interest of space the observations 
are not reported here. Next, the advantages offered by the 
electrical structure-based PMU placement over the topology- 
based approaches is demonstrated via state estimation and fault 
detection algorithms. 

IV. Applications using PMU measurements 
A. Static state estimation 

In this section, results of incorporating PMU measurements 
into a static state estimation algorithm is presented. The mini¬ 
mal number of PMUs and their optimal locations are decided 
by the solution to the PMU placement problem 0. The 
implications of both topology-based and electrical structure- 
based approaches to PMU placement on the performance of 
state estimation are discussed. In static state estimation, the 
control center accumulates a single snapshot of measurements, 
collected from meter readings and other available informa¬ 
tion from supervisory control and data acquisition systems 
(SCADA), to provide estimates of the system state. The 
measurement model is first presented, into which the PMU 
measurements are augmented. 












For a power grid comprising N buses, the system state vec¬ 
tor to be estimated is s = [02, ■ ■ ■ ,0iy,Vi,..., Vn]- The con¬ 
trol center accumulates a measurement vector y € > 

2N — 1 to estimate the state vector s. The measurements 
consist of active and reactive power flows and power injections 
at the buses. The active and reactive power flows between 
buses i and j are denoted Pij and Qij, respectively, while 
Pi and Qi denote the active and reactive power injections at 
bus i, respectively. The following relations characterize the 
measurement function mil: 


p 

JTij 

= cos{9i - 9j + aj) 



cos(ay ) -f V^jsj cosiasj), 

(9) 

Qij 

= yiyjl%j sin(6*^ - Qj + aij) 

-^%lij + Vi^^j sin(as_^-), 

B 

— ^ ^ COs(02 Oj 

(10) 

Pr 

(11) 

Q^ 

B 

(12) 




where 7 ij and are the magnitude and phase angle of the 
admittance from bus i to bus j, respectively, and jsj and asj 
are the magnitude and phase angle of the shunt admittance 
at bus i, respectively. = 0 if the buses i and j are 
not connected, which models the topology constraints of the 
network. The vlaues 'jij, ^sj and asj are assumed known. 

The measurement model is given hy y = f{s) + v, where 
f denotes the set of L nonlinear functions specified by (|9]l 
- (fTSli relating the state and measurement variables, and v 
denotes the L x 1 measurement noise vector assumed to be 
Gaussian with covariance matrix Y. The Jacobian matrix J of 
the network is given by the partial derivative J = . The 

weighted least-squares (WLS) estimate of the state vector s at 
the it + 1)* iteration is Jj+i = s* -f X* {y - f{st)) 

(see||22l), where the gain matrix is Lt = Jt) 

When [y ~ f{st)) attains a predefined value, the 

estimator is said to have converged, and the error covariance 
matrix of the estimate s is given by Rg = j) 

In II 23 I . a procedure was devised to incorporate the PMU 
measurements into the state estimation framework described 
above. By using the fact that the PMU is designed to record not 
only the voltage magnitude and phase angles at the bus where 
it is installed, but also current phasors of all lines incident on 
that bus, the measurements at the PMUs is defined by y = 

[Treat Tjmg Xeal -^img]? where Vreah Tjmg, Xeal and limg the 

real and imaginary parts of the voltage and current phasors, 
respectively, given by 


^,real 

= [Uilcos^i, 

(13) 

^,img 

= |U|sin6>i, 

(14) 

-^2 j, real 

= (|Ui| cos9i - \Vj\cos9j)gij - 

(|17,|sin6>,) 


- (|Ui| sin^i - \Vj\sm9j) bj, 

(15) 

hjjxng 

= (1 Vi 1 cos - 1 Vj 1 cos 9j ) bij P hjo 

(|Ui| sin 9i) 


-f (ll/ilsin^i - \Vj\sin 9j) gij, 

(16) 


where and bij are the conductance and susceptance be¬ 
tween buses i and j, respectively. The measurement model for 
the PMU measurements is, therefore, given by y = f{s) -f v, 
where / is defined by (fOT l - ([TbT l. and v is the PMU 
measurements noise vector assumed to be i.i.d. Gaussian with 
covariance matrix Y. The new measurement model incorpo¬ 
rating the PMU measurements is given by 


y 

y 


/(s) 

/(s) 


(17) 


The corresponding Jacobian matrix is given by 


d(f(s)) 1 

ajs) 

oifM) 

dis) 


(18) 


The error covariance matrix of the new estimate Snew is 
= ('^n^w'*'rrew>^new) , where the noise covariance 


matrix 


Y 0 
0 Y 


The WLS estimate of s at the (f -|- 1)* iteration, with PMU 
measurements, is Sj+i = St 


T fT "Y— 1 


Xnew.t (y - fist)) + 

(y - where the gain matrix inew.t = 


J/Y-iJt +J/Y-iJt 


, and the Jacobian matrix J is 


given by the partial derivative J = ■ 

B. Fault detection 


As before, we model a power system with N buses and 
K branches modeled using an undirected graph Q = (V,i£), 
where V = W} denotes the set of buses and £ = 

{(it, ji),..., (tif, jif)} denotes the set of transmission lines. 
We assume that no bus is connected to itself. As an example, 
we consider the IEEE-14 bus network with 14 buses and 21 
branches. The voltage phasor angles at time t are collected 
in the A^-length vector x{t), t = 1,..., T. The active power 
flow between the buses at time t is given by the N x N skew 
symmetric matrix Y{t), t = 1,... ,T. We denote by B the 
N X N susceptance matrix on the transmission lines; Bij = 
Bj j is the susceptance between the buses i and j if {i,j) € £ 
and Bij = 0 otherwise. We employ the linear DC power flow 
to describe the energy flow through each transmission line; 


Yij{t) = Bij [xi{t) - Xj{t)] ,t = l,...,T. (19) 

We let the iC-length vector b with elements Bij for 
ii,j) € £ and i > j. Similarly, the elements Yij{t) for 
{i,j) C £ and * > j are collected in the iC-length vector 
zif). Thus, ([T 9 I 1 can be expressed as z{t) = diag{5}I9a:(f), 
where diag{5} is the KxK diagonal matrix with the elements 
of the vector b along its main diagonal. The K x N matrix D 
is defined as follows: Eor the fc* transmission line iik,jk) G £ 
and ik > jk, we let Dkj^. = 1 and Dkj^. = — 1 , while the 
other elements in the /c* row of D assigned zero. However, in 
practice, measurements of the power and the voltage phasor 
angles are noisy especially when one deals with real-time 
monitoring of the electric grid, i.e., we observe 


z{t) = zit) p w^{t), 
x{t) = x{t) + w^{t), 


( 20 ) 

( 21 ) 


















where and {wx{t)Y^^^ are independent Gaus¬ 

sian noise processes with Wz{t) ~ JV {0,a'^l), Wx{t) ^ 
M (O, cr^J), where 0 is the vector of all zeros, I denotes the 
identity matrix, and cr^ and cr^ are the noise variances. 

Taking observations over T time units, with a sampling 
rate of one sample/unit-time, the observation model can be 
compactly represented in matrix form as follows: 

Z = diag{b}DX, 

Z = Z + W„ (22) 

X = x + w,, 


where the dimensions of the matrices Z and are K x 
T, while those of matrices X and Wx are TV x T. This is 
commonly referred to as the linear errors-in-variables (EIV) 
model ll24l . We consider the problem of detecting changes 
in the susceptance vector b based on the noisy observations 
Z and X. Specifically, we assume knowledge of a vector bg 
corresponding to the nominal behavior of the grid and pose a 
binary hypotheses testing problem: 


Ho : b = bo, 
Hi : b ^ bo. 


(23) 


Under both hypotheses, Z and X are unknown and have 
to be estimated. A statistical model for (l22l i was proposed 
in II25II - IIZTII by letting x{t) be random variables centered 
around x{t). We obtain the maximum likelihood estimates 
of Z and X and use it in the context of detection to derive 
a hypothesis test termed TML-GLRT. The full details of 
TML-GLRT were first developed in ll28l . In the TML model, 
the observed x{t) are deterministic, while x{t) are random, 
i.e., p{x{t)\x{t)) ~ M {x{t),a^l). Conditioned on x{t), 
we have p{x(t)\x(t)-,b) ^ A/" (diag{6}£)a;(f), . 

Therefore, the marginal distribution of z{t) is 
p{z{t);x{t),b) ~ M {diag{b}Dx{t), H{b)), where 

H{b) = all + aldiag{b]D'^Dd\ag{b}. The TML-GLRT is 
given by 


txML — log ■ 


uia.yibYY^p{z{t)]x{t)\b) Hi 




^ 7tml, 


(24) 


where 7tini is a fixed threshold. txML can be written in the 
following simplified form: 


T T 

tTML=niax^logp(i(f);£(f);b) - ^ logp (£(f); i(f), bg) 

= ^-Tv{A^{bo)H-\bo)A{bo)} + +1 log |Tf()| 

-minTr{AT(b)Jf-i(b)A(b)} +^log|/T()|, (25) 

where A{b) = Z — d\ag{b}DX. The test (l25T l compares 
the test statistic hmi to the threshold 7tnii, respectively. We 
choose the hypothesis Ho if the test statistic is less than the 
threshold and Hi if it is greater than the threshold. Given 
a large number of samples, the asymptotic performance of 
TML-GLRT under the hypothesis Hq is given by 2ttini 
where Xk is the Chi-squared distribution with K degrees of 
freedom ll29l . This result is independent of the specific value 


of the unknown X. A reasonable approach to choosing the 

threshold for a given false alarm rate a is 7tmi.a = 

where is the cumulative distribution function of the 

\k 

Chi-squared distribution with K degrees of freedom. 


V. Experiments and discussion 

Eor state estimation, we consider IEEE 14, 30 and 57 bus 
systems, and 50 time-sampled intervals with successful load 
flows under different loading conditions. The measurements 
of voltage magnitudes, power injections and power transfers 
are corrupted with additive white Gaussian noise having zero 
mean and standard deviations of 0.1% and 2% of actual values 
for voltages and powers, respectively. The generator outputs 
are adjusted to avoid overloading the swing bus. 

The following three experimental scenarios are considered: 
(1) Normal load operation: the variation in loads are made to 
follow a linear trend of 10%, 20% and 30%, with a random 
fluctuation of 3% over the 50-sample time interval, without 
considering bad data in the noisy measurements accumulated 
at the control center; (2) Presence of bad data: at time instant 
t = 25, one active and one reactive power transfer element 
with random error is introduced, without any errors at other 
time instants; and (3) Sudden load change: a portion of the 
load at one PQ—bus is disconnected at time instant f = 30 to 
simulate a sudden load change scenario. We performed M = 
100 Monte Carlo simulations, and for each simulation a new 
set of noisy measurements is included. 

Eollowing the setup in ll23ll . PMU measurements are as¬ 
sumed to be obtained from devices with calibrations to correct 
for current and transformer turns ratio errors. Phasor magni¬ 
tude errors are assumed to have a standard deviation of 3% of 
the actual measurement, while phase angles are assumed to be 
in error only due to time synchronization errors. Phase angle 
error is assumed to have a standard deviation of 0.02 degree. 

We consider the following performance indices commonly 
used in power systems lIMIl : (i) The estimation error at every 
time step, averaged over the set of simulations is given by 
G = where St^m is the estimated 

vector at the to* simulation, (ii) The performance index for 
the measurement vector at time step t: 


E h L", _ .j.noisi 

m=l yt,m 


Pi 


index,t — 


noisefree I 


- 77 “ 

Z^m—1 \ yt,m Ut, 


noisefree I 


(26) 


where, at the to* simulation, yt,m is the estimated measure¬ 
ment vector, y" is the noise-free measurement vector, and 
yt,m is the noisy measurement vector at the control center. 

Two different scenarios in the experimental setup are stud¬ 
ied. In the first scenario (labeled scenario 1), the same number 
of PMUs are employed for both the topology- and electrical 
structure-based approaches. As observed in Section |III1 the 
optimal location of these PMUs are different. For instance, let 
us consider the case of IEEE 14 bus system as an example. As 
shown in Table |I] the topology-based PMU placement solution 
yields a minimum of 4 PMUs, while the electrical structure- 
based solution suggests a minimum of 7 PMUs. Therefore, 
for scenario 1, phasor measurements from only 4 PMUs are 
incorporated into the WLS algorithm for both the topology- 
and electrical structure-based PMU placement schemes. 




IEEE bus system 

WLS 

WLS with PMUs (topology) 

WLS with PMUs ES (scenario 1) 

WLS with PMUs ES (scenario 2) 

14 

5.6511 

5.6615 

4.5142 

3.6457 

30 

3.1009 

2.2388 

1.9734 

1.6535 

57 

4.7611 

3.5512 

2.2149 

1.0091 


TABLE II: Estimation error e x 10^ p.u. of WES and WES with PMU measurements, during normal operating condition. ES: 
electrical structure-based. 


IEEE bus system 

WLS 

WLS with PMUs (topology) 

WLS with PMUs ES (scenario 1) 

WLS with PMUs ES (scenario 1) 

14 

0.4402 

0.3803 

0.3391 

0.3102 

30 

0.4797 

0.4198 

0.3900 

0.3596 

57 

0.5002 

0.4293 

0.4000 

0.3201 


TABLE III: Time averaged performance index Pindex (p u.) of WES and WES with PMU measurements, during normal operating 
condition. ES: electrical structure-based. 


Note 5.1: In order to conform to the constraint of com¬ 
plete network observability on the optimization setup (|7]i, the 
number of PMUs installed on the grid is given by Table H] 
however, one only utilizes phasor measurements from 4 PMUs. 
The phasor measurements from the remaining 3 PMUs is 
neglected/discarded by the WES algorithm for scenario 1 
though they remain installed on the grid. 

Eor the second scenario (labeled scenario 2), the number of 
PMUs follows from the second and third columns of Table U 
for the topology- and electrical structure-based approaches, 
respectively. And, WES is implemented with phasor measure¬ 
ments from these PMUs. 

In Table HI] the estimation error averaged over time for the 
14, 30 and 57 bus networks under normal operating conditions 
are tabulated. As shown, the performance of WES improves 
by incorporating PMU measurements for all the bus networks. 
Of particular importance is the improvement in performance of 
WES aided by the electrical structure-based PMU placement 
(scenario 1) compared to WES incorporated with phasor data 
from topology-based placement scheme; compare values in 
columns 3 and 4. The improvement is due to better charac¬ 
terization of the connectedness between buses in the network 
as provided by the electrical structure-based PMU placement. 
The best performance of WES is observed when it is aided 
by phasor measurements from all the PMUs obtained from 
electrical structure-based placement scheme (scenario 2); this 
is tabulated in column 5. Next, in Table HIH the performance 
index given by (l26l l averaged over time for all the three bus 
systems, under normal operating conditions are tabulated. Here 
again, an improvement is seen in the performance of WES 
by incorporating PMU measurements. Scenario 1 yields a 
better performance than the topology-based scheme, however, 
scenario 2 offers the best performance. 

In Pig. [T] the performance index Pindex,* given by (l26l l is 
plotted as a function of time for 100 time units while the 
system is operating under normal load conditions. In Pig.|2| the 
plot of Pindex,t is shown when the system has bad data injected 
at time instant t = 50, while in Pig. |3| Pindex.t is plotted when 
the system suffers from sudden load change at t = 75. Eor 
all the three experiments, both scenario 1 and scenario 2 are 
considered. It is noticed that the performance of WES shows 
an improvement when phasor data is incorporated into the 
algorithm. Also, in all the case studies, the electrical structure- 
based PMU placement scheme yields a superior performance 


compared to the topology-based approach. 

To analyze the performance of TML-GLRT with PMU 
measurements, we consider the IEEE 14 test bus system. The 
noise variances are cr^ = cr^ = 0.01. We perform 10^ Monte 
Carlo simulations. At each realization, the elements of X are 
generated independently as standard normal random variables. 
The active power flow Y between the buses were obtained 
using the power flow algorithm in MATPOWER EOll : the 
susceptance vector b was then computed using Y. Under Hq, 
we use the susceptance parameters given by the test profile. 
Under Hi, we apply a change of -2% to every element of 
b. Here, -2% essentially means that every element of b is 
made smaller by a factor of 2%; this number can be arbitrarily 
chosen without loss of generality. We compute the variation of 
probability of detection (P^) for a change in the susceptance 
parameter versus time, using a fixed value of false alarm rate 
(a) uniformly picked between [0, 0.2] for which the threshold 
7tmi,a is calculated. This procedure is repeated for different 
values of a S [0, 0.2] resulting in a set of P^s versus time. 
We plot the average (over the number of a) Pd versus time. 

Eor the topology-based approach, we obtain an optimum 4 
PMUs to be installed on buses 2, 6, 7 and 9. Eor the electrical 
structure-based approach, we need an optimum number of 7 
PMUs to be installed on buses 1, 3, 8, 11, 12, 13 and 14. The 
plots of Pd versus time for one frame lasting 20 time units 
when the PMUs are placed employing both the approaches 
are shown in Pig. |4a| Eor a fairer comparison, we modify the 
experimental setup as follows: We employ the same number 
of PMUs for both the topology and electrical structure-based 
approaches to analyze the detection performance. Specifically, 
we employ only 4 PMUs transmit their phasor data to the 
control center for both PMU placement approaches for fault 
detection. More precisely, for the topology-based approach we 
allow PMUs on buses numbered 2, 6, 7 and 9 to transmit as 
before. Eor the electrical structure-based approach, we now 
allow only 4 PMUs on buses numbered 8, 11, 12 and 14 
transmit their phasor data. This choice was based on the 
electrical centrality measure between various buses on which 
the PMUs are placed; this was characterized in an earlier work 
by the authors ImI . As seen in Pig.|4bl the electrical structure- 
based approach has a higher detection performance compared 
to it topological counterpart. 

Einally, we conduct an experiment in which we allow lesser 
number of PMUs for the electrical structure-based approach 
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Fig. 1: Performance index under normal load conditions. 





(a) IEEE-14 bus system 


(b) IEEE-30 bus system 

Fig. 2: Performance index under bad data conditions. 


(c) IEEE-57 bus system 





(a) IEEE-14 bus system (b) IEEE-30 bus system (c) IEEE-57 bus system 

Fig. 3: Performance index under sudden load change. 


to transmit compared to the topology-based approach. Specifi¬ 
cally, we allow only 3 PMUS (on buses 14, 8 and 12, based on 
the based on the electrical centrality measure) to transmit for 
the electrical structure-based approach, while continuing with 
4 PMUs for the topology-based approach. The resulting plots 
of Pd versus time are shown in Fig. |43 where it is clearly seen 
that Pd for the electrical structure-based approach is slightly 
higher than the topology-based approach. 

Note 5.2: For Fig. |4b] and Fig. @3 it should be noted that 
we have installed the optimum number of 7 PMUs for the 
electrical structure-based approach, so as to achieve complete 
network observability. However, we allowed only 4 PMUs to 
transmit their phasor data to the control center. Though the 
remaining 3 PMUs record phasor data across the lines on 


which they are installed, they do not transmit them to the 
control unit. So long as phasor data from the “high priority” 
{i.e., higher electrical centrality measure) buses are obtained, 
the test statistic can be computed for hypothesis testing. 

The electrical structure-based approach provides a stronger 
characterization of the electrical influence between network 
components compared to the topology-based approach. There¬ 
fore, the probability of detection of change in the susceptance 
parameter is higher when the PMUs are installed based on the 
electrical structure-based approach compared to the case when 
they are installed using the topology-based approach. 

VI. Conclusion 

The PMU placement problem was addressed from the 
perspective of the electrical structure of the grid. PMU mea- 
















































































(a) (b) (c) 

Fig. 4: Probability of detection versus time, for one frame lasting 20 time units, for IEEE-14 bus system. 


surements were incorporated into WES used for static state 
estimation and TML-GLRT used to detect changes in the 
susceptance parameters of the grid. The main result was that 
on the one hand, a larger number of PMUs were required to 
achieve complete network observability, while on the other 
hand, it improved the performance of practical applications 
such as state estimation and fault detection. Given that the 
electrical structure provides a more comprehensive description 
of interconnections in the power network, our results promote 
a realistic perspective of PMU placement in the grid. 
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